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Abstract. Growth curves of the near-infrared (NIR) 
magnitude as a function of the aperture have been built 
and used to derive NIR total magnitudes from aper- 
ture data taken from the literature. By cross-correlating 
with optical and redshift data, absolute magnitudes and 
optical-to-NIR colors have been computed for some 1000 
galaxies of different types. Significant color gradients are 
observed, underlining that small aperture colors may lead 
to a biased picture of the stellar populations of galaxies. 

A statistical analysis using various estimators taking 
into account the intrinsic scatter has been performed to 
establish relations between the colors, the morphological 
type, the inclination or the shape, and the intrinsic lumi- 
nosity. 

The combination of the optical and the NIR should 
obviously improve our understanding of the evolution of 
galaxies. Despite the intrinsic scatter, especially among 
star-forming galaxies, optical-to-NIR colors show a very 
well defined sequence with type, blueing by 1.3 mag 
from ellipticals to irregulars. The colors of spiral galaxies 
strongly redden with increasing inclination and put new 
constraints on the modeling of the extinction. No such ef- 
fect is observed for lenticular galaxies. We also find that 
rounder ellipticals tend to be redder. 

A color-absolute magnitude relation is observed inside 
each type, with a slope significantly steeper for early and 
intermediate spirals than for ellipticals or late spirals. This 
stresses the importance of considering both the mass and 
the type to study the star formation history of galaxies. 

Key words: Galaxies: evolution - Galaxies: fundamental 
parameters - Galaxies: photometry - Galaxies: statistics 
- Infrared: galaxies ~ dust, extinction 
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1. Introduction 

Although color-magnitude diagrams have now been ob- 
tained for a few very close galaxies - most of them dwarf 
-, allowing to study directly their stellar populations (e.g. 
Schulte-Ladbeck & Hopp 1998; Aparicio 1998), the syn- 
thesis of the spectral energy distribution (SED) of nearby 
galaxies remains the most efficient and systematic method 
to trace their star formation history up to z = 0. 

Early optical studies have shown a clear correlation 
between the SED - and hence the star formation history 
- and the galaxy Hubble type. However, the age-star for- 
mation timescale and age-metallicity degeneracies make it 
highly desirable to extend the wavelength range to the ul- 
traviolet (UV) and the near-infrared. Whereas the ultravi- 
olet is related to the young stellar populations in currently 
star- forming galaxies, the near-infrared is dominated by 
old giant stars and is a measure of the star formation rate 
integrated from the beginning, i.e. of the stellar mass of 
the galaxy. When compared to shorter wavelengths, the 
NIR may moreover put constraints on the stellar metal- 
licity, owing to the high sensitivity of the effective temper- 
ature of red giants to the abundance of heavy elements. 
In combination with optical data, it finally provides clues 
on the amount and the distribution of dust because of the 
different extinction at both wavelengths. 

Observations in the UV and the NIR are however ham- 
pered by the atmosphere. Though interesting to study 
the detailed star formation history of a specific galaxy, 
the spectra obtained from Space are too scarce to afford 
a statistical analysis from which the general trends as a 
function of type or mass, or the effect of the dust, may be 
derived. In spite of the loss of spectral resolution, broad- 
band colors in the J, H and K atmospheric windows are 
more promising. As the spectra, they are usually obtained 
in small apertures not representative of the whole galaxy 
because of the color gradients, especially for disk galaxies, 
and have to be extrapolated to be compared to the op- 
tical data. Such extrapolation is however much easier for 
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colors than for spectra and can be applied to hundreds of 
objects. 

In the following, we analyze a catalog of infrared aper- 
ture magnitudes (section ^ and combine it with optical 
catalogs to derive NIR growth curves of the magnitude as 
a function of the aperture (section ||) . We then compute 
total and effective magnitudes and colors (section The 
NIR and optical-to-NIR colors are analyzed statistically as 
a function of type, luminosity and inclination in section ^j. 
We finally discuss our results in section ^. 



type 


fee 


kj 


fen 




E 


5.1z 


-0.3z 


-0.l2 


-2.82 


SO 


4.92 


-0.3z 


-O.I2 


-2.82 


Sa 


4.O2 


-0.32 


-O.I2 


-2.82 


Sb 


3.62 


-0.32 


-O.I2 


-2.82 


Sc 


2.62 


-0.3z 


-O.I2 


-2.82 


Sd 


2.2z 


-0.52 


-0.32 


-2.82 


Im 


l.8z 


-0.72 


-0.42 


-2.82 



Table 1. fc-corrections in B, J, H and X as a function of 
the type. 



2. The optical and infrared data 



The near-infrared data use d in t his paper come from the 
CIO catalog (Gezari et al. 1997 ) which is a compilation 
of the NIR observations published before 1995. They are 
given in a large variety of systems (magnitudes toa, 
F^, etc.) and units, and are first converted to magnitudes 
assuming Vega has a magnitude of 0.03. The data are 
reduced from the wavelength A to the central wavelengths 
of the J (Aj = 1.24^m), H (Ah = 1.65££m) and K (Ak = 
2.21/im) filters of Bessel & Brett ( 1988 ), assuming the 
following color equations: 



m,^ = mx + (Aj — A) 



<J-H> 
A.T — Ah 



mn = mx + (Ah - A) 



<J-H> 
Aj — Ah 



if A < 1.65/im, 



iTT-x + (Ah — A) — ^ ^ — — if A > 1.65/j,m, 



niK = mx + (Ak - A) 



Ah — Ak 

< H -K > 
Ah — Ak 



where < J - i? 0.7 and < H - K >= 0.2 are the 
typical J — H and H — K colors of a normal galaxy and 
depend very little on the galaxy type or the aperture. 

The CIO catalog is cro ss-co rrelated with the RC3 cata- 
log (de Vaucouleurs et al. 1991 ) to correct for the Galactic 
extinction. The NIR extinction is computed from Ab as 
follows: Aj = 0.182^B, = 0.098^b, Ak = 0.068Ab. 

The magnitudes are also corrected for the redshift z 
as mx{0) = mxiz) — kx{z). Values of kx{z) are given in 
Table |^ for different types of galaxies and have been com- 
puted using the PEGA SE m odel of spectral evolution (Fioc 
& Rocca-Volmerange 1997 ). They should be used only at 
z < 0.1. The redshift is taken from the NED database. 
For a few galaxies, z is unknown and we compute the 
fc-correction (but not the absolute magnitudes) assuming 
z = 0.01, the mean redshift of our catalog. 

The RC3 provides us also with the morphological 
type T, the ratio R25 of the major axis {D25) to the 
minor axis of the ellipsis corresponding to the isophot 
fi = 25mag/arcsec^, the total i?-magnitude Bt extrap- 
olated to an infinite radius and the circular effective aper- 
ture Ae containing half the light emitted in the _B-band. 



When available, we prefer to take Bt and Ac from 
the m ore recent catalog Hypercat (Prugniel & Heraudeau 
1998). This catalog also gives the photometric type Tp 



corresponding to the best-fitting growth curve B{X,Tpj 
of the B-magnitude as a function of the circular aperture 
A, where X = \ogiQ{A/Ac) (see Appendix A). 

3. The NIR growth curve 

The detailed 2D-fitting of the profile of galaxies (de Jong 
& van der Kruit 1994) is certainly the best way to com- 
pute their total magnitudes, but surface photometry in the 
NIR is available for too few of them to perform a statis- 
tical analysis. The largest such stu dy is b a sed on only 8 6 
spiral galaxies of all types (de Jong 1996a , 1996b , 1996c ), 
to compare to typically 100 galaxies per type in this pa- 
per. We therefore prefer to extrapolate aperture magni- 
tudes with a growth curve. The computation of total (i.e. 
asymptotic) NIR magnitudes comparable to the optical 
ones determined by, e. g., de Vaucouleurs et al. ( 1991 ) or 
Prugniel & Heraudeau ( 1998| ) is however made difficult by 
the small apertures achieved at these wavelengths. Few 
galaxies have enough data from small to large apertures 
to constrain the shape of the growth curve and one has to 
use another method, combining observations of different 
galaxies having presumably the same curve (cf. Griersmith 
1980). To this purpose, one needs to scale the apertures to 
a characteristic length Ac for each galaxy. The NIR growth 
curve A4{A/Ac) is then built by plotting ■m{A) ~ m{Ac) as 
a function of login (A/A c)- Usually (Gavazzi et al. |1996a ; 
Gava zzi fc Boselli 1996| ; Tormen & Burstein 199£| ; Frogel 



et al. 1978; Aaronson et al. 197S| ), the D25 diameter has 



been used as characteri stic le ngth. As noticed however by 
de Vaucouleurs et al. ( 1976 ), the D25 diameter depends 
not only on the shape of the profile but also on the cen- 
tral surface brightness. For this reason, the ratio of the 
effective aperture Ae to D25 is not a constant. 

We prefer to adopt here an other procedure. We may 
expect some relation between the optical profile and the 
NIR profile. For this reason, we decide to build a growth 
curve as a function of the B photometric type Tp taken 
from Hypercat and of X = \ogiQ{A/ Ac). The infrared 
magnitude m{Ac) in the B effective aperture is computed 
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by interpolation or by a small extrapolation in log]^Q(A). 
Rather than a straight line, we have used the functions 
proposed to compute H magnitudes at yl = D25 by 



Gavazzi et al. (1996a): 

M{X25) = a + fe(-0.35X25 + I.68XI5 



0.07X1,) 



where X25 — \ogiQ{A/ D25), and fitted a and b to the 
observations. This takes into account the curvature of the 
growth curve near and improves the determination of 
m{Ac). Because they are polynomials, these functions are 
however not suitable to extrapolate the magnitude to the 
infinity. 

The J, H and K data have been combined to build the 
growth curve. This is justified by the fact that only small 
J — H and H — K color gradients are observed (Aaronson 



1977| ; Frogel et al. |197^ ). The growth curves are plotted 
in Fig. |l| for 7 bins of photometric type. As evidenced 
by this graph, the NIR growth curve is flatter than the 
optical one, especially for Tp corresponding to interme- 
diate spirals: the NIR emission of galaxies is more con- 
centrated than in the optical. This behavior is expected 
because the bulge is redder than the disk and is there- 
fore more prominent in the NIR than in the optical. The 
scatter is also higher for intermediate types, which might 
be due to an inclination dependency of the growth curve 
(Christensen |1990 ; Kodaira et al. 1990). We have therefore 
distinguished in Fig. |l| between galaxies with i?25 < 0.3 
(rather face-on for disk galaxies) and i?25 > 0.3, but no 
obvious trend is observed and we will neglect any inclina- 
tion dependency in the following. 

Because of the scatter and the lack of data at high 
aperture, where the flattening of the growth curve puts 
constraints on the curvature, attempts to fit growth curves 
similar to those of Prugniel & Heraudeau ( 199§| ) with 3 
parameters (mx, \ogn){Ac) and Tp) did not prove success- 
ful. Plotting the NIR magnitudes versus the B magnitudes 
however reveals a nearly linear - though scattered - rela- 
tion, between these quantities, suggesting to adopt a NIR 
growth curve of this form: 

M{X,Tp)^sB{X,T^). 

This ensures that the extrapolation at infinite A converges 
since limx^oo S(X, Tp) = 0. Such relation has been fitted 
for each photometric type bin (Fig. ^ dashed lines) at 
X > -1.5. 

The dispersion is due both to the uncertainties in 
the individual data and to the intrinsic scatter in the 
shape of the NIR growth curves for a given B photomet- 
ric type. The following convention is adopted hereafter: 
X M{^,,cf'^) means that x is distributed according to 
( or the density probability of x is) a Gaussian with mean 
H and variance . Let us assume that, for each galaxy, 
the distribution of individual data around the best-fitting 
NIR growth curve sB{X, Tp) is - Af{0, cr^ J and that the 
distribution of the parameter s characterizing the growth 



curve is s ~ J\f{so, cr^^). At any aperture, 

m{X) - m(0) ~ M{so [BiX, Tp) - B{0, Tp)] , 

al^+al[B{X,T,)-B{0,T^f}. 

A maximum likelihood estimation of these parameters 



yields 



0.06, nearly independent of Tp, — 0.10 



for early types (Tp < —1) and cts,, ~ 0.15 at later types. 
The value of sq we obtained is plotted as a function of the 
median Tp of each bin in Fig. ^ and may be approximated 
with the following formulae: 



So 



0.87 
0.74 
0.64- 



if Tp < 
0.029Tp 
O.OllTo 



-4.5, 
if 



4.5 < Tp < 2.5, 



if Tp > 2.5. 



The mean value sq is less than one for all types, corre- 
sponding as expected to a blue-outwards gradient. This 
gradient is low for early- types, increases for spirals, peak- 
ing at Tp G [2 — 4] (~ Sb) in good agreement with what has 
been obtained at optical wavelengths by de Vaucouleurs & 
Corwin (1977) and remains constant or slightly decreases 
for late types. The typical blueing of the B — NIR color 
of spirals from the effective aperture to the infinity is 
—0.75(1 — s) ~ —0.2, close to the median value (—0.19) 
in B — H and B — K we have computed from the profiles 
published by de Jong ( |1996a[ ). 

4. Effective and total near-infrared magnitudes 

J^.l. Apparent magnitudes and colors 

Effective and total near-infrared magnitudes have been 
computed for all the galaxies with known effective aper- 
ture by fitting the growth curves sB{X^ Tp) introduced in 
the previous section to the observed data {Xi, mi,ai), i € 
[l,n], where ai is the uncertainty on the magnitude 
rui. The asymptotic magnitude depends primarily on the 
shape of the growth curve at large aperture. To avoid 
problems with small apertures suffering from seeing or off- 
centering, or which are contaminated by nuclear emission 
or a central starburst, we remove all the observations at 
X < -1 from the fit. 

The procedure to compute the total and effective NIR 
magnitudes mx, and the corresponding uncertainties 
(T(mT), <y{mc) due to the uncertainties in the aperture 
magnitudes is detailed in Appendix B. 

The global uncertainty on mx taking also into 
account the uncertainties in Tp and A^ is finally computed 
from 



and the same for . 



a^{A,) 



4 



M. Fioc & B. Rocca-Volmcrangc: Optical-to-NIR colors of nearby galaxies 







-1.5 -1 -0.5 0.5 1 -10 -5 5 10 



i°gio(VAe) Tp 

Fig. 1. (a) to (g): Growth curve of the NIR magnitude as a function of the aperture for the different photometric 

types. Hollow and filled triangles correspond respectively to galaxies with logj^g R25 < 0.3 and logj^Q R25 > 0.3. The 
solid and dashed lines are the mean optical and NIR growth curves. They are normalized to the effective aperture in 
the B-band and to the corresponding magnitude, (h): The circles and vertical segments are the mean value so and 
the intrinsic scatter (T.,,, of s as a function of the median value of Tp in each bin. The broken line is the 3-slope fit to 
so(Tp) given in the text. 
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In many cases, the photometric type is unknown and 
we estimate it from the morphological type with the fol- 
lowing relations derived by cross-correlating Hypercat and 
the RC3: 



Tp = -2.3-(-0.5T 



if T < -1, 
= -1.32 + 1A8T if - 1< T < 4, 
= 1.4 + 0.8r if T > 4. 

We then simply assume cr{Tp) — 3. 

The uncertainty on the colors mi 
(mi,TO2) G [J,H,K]'^, is computed as 



m2, 



where 



cr^(TOi) + cr^(m2) -f 



9(toi — 7712) 



d{mi 



m.2) 



However, this is clearly an overestimate because J, H and 
K data have usually been observed simultaneously and 
their errors are presumably highly positively correlated. 

When B is one of the band, e.g. mi = Bt, the partial 
derivatives are unknown. We then estimate <^rni-m2 ^-s 



This will usually be a slight overestimate for the partial 
derivatives of Bt and m2 should have the same sign, but 
most of the uncertainty comes actually from the NIR data 
themselves. 

4-2. Absolute magnitudes 

The absolute magnitude Mx is computed as 

Mt = niT - 25 - 5 logiQ(t;v/i?o) 

where vy is the velocity of a galaxy derived from the red- 
shift z and corrected for the movement of the Sun in the 
restframe of the Virgo cluster according to the eq uation s 
(18), (19), (20), (31) and (32) from Paturel et al. ( |1997D . 
We assume Hq = 65km.s~^.Mpc~^. The uncertainty on 
the absolute magnitude is computed from 

' {dvv/dz)^a^z)+vl 



- 



In(lO) 



where Vp = 350 km. s is the typical peculiar velocity 
of galaxies in the Las Campanas Redshift Survey (Lin et 
al 



1996) 



At low redshift, peculiar velocities perturb the 
redshift-distance relation and increase the uncertainty on 
the absolute magnitude. All the absolute magnitudes of 
the galaxies with [{dvy /dz)^a'^{z) + v^]^^^ > vy/S have 
been discarded in the following. <7{z) has been computed 
according to the formulae (1) and (2) of the Appendix 
A of Paturel et al. ( |1997[ ), taking into account both the 
internal uncertainty on the values of the redshift and the 
discrepancies between them. 

Covariances between the colors and the absolute mag- 
nitudes have also been computed. 



type 


T (RC3) 


E 


r < -3.5 


SO 


-3.5 < T < -0.5 


Sa 


-0.5 < T < L5 


Sb 


1.5 < T < 3.5 


Sbc 


3.5 < T < 4.5 


Sc 


4.5 < T < 5.5 


Sd 


5.5 < T < 8.5 


Im 


8.5 < T 



Table 2. Correspondence between the types used in this 
paper and the RC3 types. 



5. Statistics of optical-near infrared colors 

5.1. J — H , H — K and J — K colors 

The galaxies have been distributed in eight broad types 
according to the morphological type T taken from the RC3 
(see Table |). The J ~ H, H - K and .J - K total and 
effective colors have been plotted in figure || as a function 
of the type. 

Because the uncertainties computed ior J — H, H — K 
and J — K are overestimates, we simply give in table ^ the 
median color for each type and compute its uncertainty by 
a bootstrap with replacement. The total and effective col- 
ors show almost no difference. Whereas H — K seems to 
be nearly independent of the type, J — H and J — K de- 
crease clearly at the late types, indicating that the J-band 
is contaminated by young stars. The colors are also redder 
for early-type spirals, which is difficult to understand with 
current models of stellar evolution and atmospheres. This 
effect is observed in both the total and effective colors. 
The extrapolation to compute the latter being small, this 
may not come from a problem with the growth curve. The 
extinction, negligible in the NIR, is also unlikely to cause 
this phenomenon. A possible explanation is that we begin 
to see in H and K the tail of the dust infrared emission at 
short wavelengths. We note in particular that the J and K 
samples contain a non-negligible fraction of galaxies with 
nuclear activity of which central colors are in some cases 
as red as J — H^H — K^l. 

5.2. B~J,B-HandB~K colors 

The catalog is large enough to perform a statistical anal- 
ysis of the colors as a function of the type, the luminosity 
and the inclination or the shape (-R25) of the galaxy. An 
interesting quantity is also the intrinsic scatter in the col- 
ors at a given type, luminosity and i?25, for it is a measure 
of the variations of the star formation history and of the 



effects of dust (e.g. Peletier & de Grijs 1998; Shioya & 
Bekki 1998] ) . Note however that the intrinsic scatter may 
depend on the type-binning, especially if the colors evolve 
rapidly from type to type. 
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Fig. 2. Total and effective NIR colors as a function of type. A random component in [—0.5, 0.5] has been added to the 
type to make the graph clearer. The diamonds are the median colors per type. 
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type 


iV 


iJ-H)T 


(J-H), 


iV 


{H - K)t 


{H - K), 


N 


{J-K)t 


(J - K), 






median 


median 




median 


median 




median 


median 


E 


140 


0.71 ± 0.01 


0.71 ± 0.01 


141 


0.21 ± 0.01 


0.22 ± 0.01 


140 


0.92 ± 0.01 


0.92 ± 0.01 


SO 


154 


0.72 ±0.01 


0.72 ±0.01 


157 


0.23 ±0.01 


0.23 ±0.01 


155 


0.95 ±0.01 


0.95 ±0.01 


Sa 


96 


0.74 ±0.01 


0.75 ±0.01 


96 


0.28 ±0.01 


0.28 ±0.01 


97 


1.01 ±0.01 


1.02 ±0.01 


Sb 


93 


0.78 ±0.01 


0.79 ±0.01 


95 


0.26 ±0.01 


0.27 ±0.01 


95 


1.03 ±0.01 


1.04 ±0.01 


Sbc 


46 


0.77 ±0.02 


0.77 ±0.02 


47 


0.23 ±0.02 


0.24 ± 0.02 


47 


1.01 ±0.02 


1.01 ± 0.01 


Sc 


46 


0.72 ± 0.02 


0.73 ±0.02 


46 


0.22 ±0.02 


0.22 ±0.02 


46 


0.96 ±0.02 


0.96 ± 0.02 


Sd 


26 


0.72 ± 0.04 


0.71 ±0.04 


24 


0.22 ±0.04 


0.22 ±0.04 


24 


0.87 ±0.05 


0.88 ± 0.04 


Im 


22 


0.57 ±0.05 


0.56 ±0.03 


20 


0.24 ±0.05 


0.25 ± 0.04 


25 


0.73 ±0.05 


0.74 ±0.06 



Table 3. Median total and effective J — H, H — K and J — K colors (with their lu-uncertainties) as a function of 
type. 



The B — H sample is by far the largest (more than 
900 galaxies in B — H against 600 in B ~ J and B — K). 
It is also the most accurate (the uncertainty is typically 
0.16 mag in {B - H)t and 0.09 mag in {B - H)^) and 
the most complete in type, whereas the B — J and B — K 
sample are strongly deficient in the latest types. For these 
various reasons, we will mainly focus our discussion in the 
following on the analysis of the B — H data. 

We adopt here the formalism proposed by Akritas 
( |1999D . Let us assume that the true color Yi of the i"^ 
galaxy depends linearly on p quantities Xij,] G (the 
Xij being, e.g., log]^Q(i?25) or the absolute magnitude): 



where ^.j , the median of Xij , is introduced merely to iden- 
tify /3o with the typical color of the sample, and ti is the 
deviation from the relation due to the intrinsic scatter. 
We assume that ^ A/'(0, c^) independently of the other 
parameters. Let us write /3p+i = ct to remind that the 
scatter is also to be determined. 

The observed variables, Xij and yi, are related to the 
true variables by 



X, 



1 1 



where ~ A/'(0, erf ) and eyi ~ A/'(0, crf-j). We obtain 



P 



Assuming that Ci is not correlated with the eij and ey-i, 
and that the covariances between e^j and Cik and between 
eij and Cy-y are 7yfc and 7ijF, respectively, we obtain that 

V^-Po- Ej = l PA'^rj - Mj) ~ A/'IO, <J^^ + P^+i) with 

p p p 

J=l j=l j,k=l 



We have tested two procedures to estimate the (3j from 
our sample. The first one estimates the (3j by maximiz- 
ing the logarithm of the likelihood A with respect to the 
f3j,j G [0,p± 1], where 




P+i 



ln{(7i 



according to the formulae given above. The covariance ma- 
trix of the (3j is computed by inverting the curvature ma- 
trix {uij), e [0,p+ 1]^, where 

_ 52 [- ln(A)] 



(Kendall & Stuart |l979D 



Maximum likelihood (ML) estimators are often biased and 
we therefore wish to compare to another method. 

The second procedure (MCES estimators) has been de- 
veloped by Akritas & Bershady ( |1996D and Akritas ( |1999D . 
According to the authors, it yields unbiased estimators. 
However, the cry; are not used in the expression of the 
(3i,i < P, which means that they have all the same weight. 
This is especially not satisfying if the intrinsic scatter is 
small or comparable to the uncertainties in the colors - 
which happens in particular for elliptical galaxies ~ and 
may give an excessive importance to outliers. The authors 
also do not provide the intrinsic scatter /3p+i and we esti- 
mate it from 

^ n p 

but when the intrinsic scatter is smaller or of the same 
order than the uncertainties in the observables, it is un- 
derestimatedn and we may even obtain a negative (and 



^ This is shown by the fact that the ML and MCES esti- 
mators of the scatter in the effective colors (which are less 
uncertain) are similar and are also close to the ML estimator 
for the total colors, whereas the MCES intrinsic scatter in the 
total colors is smaller. 
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meaningless) value. We then assume (3p+i = 0. Uncertain- 
ties on the Pj have been computed by a bootstrap with 
replacement on the (x^ , yi) rather than by using the cum- 
bersome formulae proposed by Akritas & Bershady (1996) 
and Akritas ( 199g| ), which anyway are not available for 
Pp+i- 

Note that the uncertainties determined either from the 
curvature matrix or by bootstrap become themselves very 
uncertain when the number of galaxies is less than about 
20 to 30. 

In the following, values determined by the ML and 
MCES estimators are written in boldface and italic char- 
acters, respectively. 

Because we are mainly interested in "normal" galax- 
ies, we have to re ject a few galaxies with unusual colors 
(Buta et al. 1994). Most of them are extremely blue and, 
according to their optical colors (when available), are pre- 
sumably starbursting galaxies. Some blue compact dwarfs 
may also be misclassified from their morphology as ellip- 
ticals. In the opposite, our sample contains some galaxies 
with very red colors, especially in the X-band and their 
NIR emission is dominated by an active nucleus. To reject 
these galaxies, we apply an iterative method. Rather than 
using the standard deviation around the best fit, which 
is very sensitive to the outliers, we estimate the observa- 
tional scatter from the median of the absolute deviations. 
At each step, we compute 



1.48 median 

ie[l,n] 



from the sample (including the "abnormal" galaxies), 
the factor 1.48 ensuring that in the case of a per- 
fectly Gaussian distribution, the standard deviation 
is S; we define "normal" galaxies as those having 



< 3(5 and then determine 



the (3j from them only by either of the methods discussed 
above. We iterate this procedure 4 times, but the conver- 
gence is usually achieved at the third one. 

5.3. The type-color relation 

The total and effective B — J, B — H and B — K colors 
have been plotted as a function of the morphological type 
on Figure |3| A very impressive trend to bluer colors with 
advancing type is observed. The mean (J5 — H)t of irreg- 
ular galaxies is 1.34 mag bluer than the one of ellipticals, 
to be compared to 0.68 inU - B, 0.47 in B -V, 0.46 in 
V — Ic and only 0.23 in V — Rc (Fioc & Rocca-Volmerange, 
in preparation). A similar gap between the {B — H)t of 
ellipticals and irregulars was obtained in {B — H)c by Buta 
( |1995D from a sample of 225 galaxies. 

The {B — H)t colors determined by the ML estimator 
are almost always systematically redder by ~ 0.02 mag 
than those computed using the MCES estimator. This bias 
comes from the fact that the uncertainties arc smaller for 



brighter galaxies in the NIR, which are also usually redder. 
The ML estimator giving more weight to the data with 
smaller uncertainties, it is biased to the red. 

A significant scatter is also observed within each type. 
Part of it is due to the observational uncertainties, but it 
comes mainly from the intrinsic scatter in the colors. The 
intrinsic scatter increases from 0.13/ .08 , depending on 
the estimator, for ellipticals to ^ 0.4 for Sb and remains 
nearly constant at later types. 

5.4- The R25- color relation 

The color-inclination relation is potentially a very power- 
ful constraint on the amount of dust and its distribution 
relatively to stars. As a disk becomes more and more in- 
clined, its optical depth increases and the colors redden. 
This must be especially striking when one of the band 
(e.g. B) is heavily extinguished whereas the dust is almost 
transparent in the other one (e.g. H). For these reasons, 
various studies have tried to determine a color-inclination 
relation, regrouping the galaxies either as a function of 
type (e.g. Boselli & Gavazzi 1994) or as a func tion of their 
NIR absolute magnitude (TuUy et al. |l998|) . The sam- 



ples used in these studies where however small (about 100 
galaxies) and it is worth to look at the relation once again 
with our much larger catalog. 

For an oblate ellipsoid, which is a standard model for 
a galaxy disk (Hubble 1926), the inclination i (face-on 
corresponds to z = 0) is related to the true ratio qo of 
the minor to the major axis and to the apparent ratio 
q ~ l/i?25 by 



cos^ i 



q- 



,2 • 



Practically, we prefer to establish a relation between 
the colors and log;^o(^25), both because i?25 is the directly 
observed quantity and because the linear estimators used 
here provide a better fit when log]^o(^25) is used rather 
than i or cos(i). 

The relations are plotted in Fig. ^ and statistical esti- 
mators of the slope and the color are given in table ^ for 
B — H . The slopes are represented as a function of type in 
figure ^. Also shown on this graph are the slopes derived 
from the larger B — V and U — B samples. 

The typical slope of the {B — i/)T-logig(i?25) relation 
is about 1 for spiral galaxies, to be compared to 0.2-0.3 
in B - V andU - B. 

Surprisingly high is the value of the slope for Sa galax- 
ies. However, the slope of the {B — H)^- inclination relation 
would be shallower for Sa since go is a decreasing function 
of type (Bottinelli et al. 1983), i.e. Sa galaxies are thicker 



than later types, as also indicated by their smaller median 
logio(-R25)- 

A dip is observed for Sbc galaxies, which seems strange 
since the maximum extinction is ex pected precisely for 
this type (Fioc & Rocca-Volmerange 1997). Though not 
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Fig. 3. Total and effective optical-to-NIR colors as a function of the type. The diamonds are the standard colors 
determined by ML. The vertical segments show the intrinsic scatter and the horizontal ones the width of the type-bin. 
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Fig. 4. {B — -ff)T-logio(-R25) relation per type. The solid and dashed lines correspond respectively to the ML and 
MCES estimators. 
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type 


iV 


A 


[{B- 


J)t] 


a 


A)[(B- 






a 


E 


140 


3 


07 ± 


0.02 





13 


± 





02 


3 


15 ± 


0.01 







13 ± 





01 






3 


± 


0.02 





10 


± 





03 


3 


13 ± 


0.02 







± 





02 


SO 


155 


3 


01 ± 


0.02 





17 


± 





02 


3 


10 ± 


0.02 







18 ± 





01 






2 


97 ± 







i5 


± 





03 




07 ± 


0.02 







24 ± 





02 


Sa 


97 


2 


74 ± 


0.04 





31 


± 





03 


2 


91 ± 


0.04 







31 ± 





03 






2 


7^ ± 


0.04 





^5 


± 





05 




93 ± 


0.03 







27 ± 





03 


Sb 


95 


1 


63 ± 


0.04 





36 


± 





04 


2 


79 ± 


0.04 







37 ± 





03 






2 


62 ± 


0.05 





33 


± 





0^ 




19 ± 


0.0^ 







55 ± 





0^ 


Sbc 


47 


1 


48 ± 


0.06 





35 


± 





05 


2 


67 ± 


0.06 







35 ± 





04 






2 


45 ± 


0.06 




i7 


± 





10 




6ff ± 


0.05 







26 ± 





07 


Sc 


46 


2 


27 ± 


0.09 





53 


± 





07 


2 


44 ± 


0.10 







61 ± 





07 






2 


24 ± 


0.0« 





45 


± 





iO 


2 


43 ± 


0.10 







56 ± 





08 


Sd 


26 


1 


96 ± 


0.09 





37 


± 





07 


2 


16 ± 


0.09 







38 ± 





07 






1 


95 ± 


0.0« 







± 





12 




± 


0.0,9 







32 ± 





10 


Im 


28 


1 


61 ± 


0.09 





39 


± 





07 


1 


74 ± 


0.10 







50 ± 





08 






1 


52 ± 


0.10 







± 









12 ± 


0.10 







41 ± 





o,s 



Table 4. [B — J)t and (B — J)o colors per type. The values computed with the ML and MCES estimators are written 
respectively in boldface and italic characters, a is the intrinsic scatter. 



type 


iV 


/3o 




H)t] 






HU 




E 


142 


3 


76 ± 


0.02 


0.13 


± 





01 


3 


86 ± 


0.02 





14 ± 





01 






3 


73 ± 


0.02 


0.08 


± 





0^ 


3 


± 


0.02 





13 ± 





02 


SO 


157 


3 


72 ± 


0.02 


0.17 


± 





02 


3 


85 ± 


0.02 





17± 





01 






3 


69 ± 


0.02 


O.i^ 


± 





03 


3 


82 ± 


0.02 





i9 ± 





02 


Sa 


110 


3 


57 ± 


0.04 


0.32 


± 





03 


3 


73 ± 


0.04 





34 ± 





03 






3 


54 ± 


0.04 


0.31 


± 





0^ 


3 


73 ± 


0.0^ 





33 ± 





03 


Sb 


171 


3 


53 ± 


0.03 


0.41 


± 





03 


3 


73 ± 


0.03 





39 ± 





02 






3 


52 ± 


0.03 


0.^0 


± 





03 


3 


11 ± 


0.03 





40 ± 





03 


Sbc 


110 


3 


33 ± 


0.04 


0.34 


± 





03 


3 


53 ± 


0.04 





35 ± 





03 






3 


30 ± 


0.03 


0.29 


± 





0^ 


3 


5i ± 


0.03 





32 ± 





03 


Sc 


125 


3 


24 ± 


0.04 


0.44 


± 





03 


3 


42 ± 


0.04 





44 ± 





03 






3 


20 ± 


0.0^ 


0.42 


± 





0^ 


3 


39 ± 


0.0^ 





43 ± 





03 


Sd 


135 


1 


94 ± 


0.04 


0.45 


± 





03 


3 


13 ± 


0.04 





46 ± 





03 






2 


93 ± 


0.04 


0.^5 


± 





03 


3 


12 ± 


0.04 





47± 





03 


Im 


42 


2 


41 ± 


0.08 


0.46 


± 





06 


2 


56 ± 


0.08 





46 ± 





06 






2 


38 ± 


0.0« 


0.42 


± 





07 


2 


5^ ± 


0.0,? 





± 





06 



Table 5. {B — H)t and {B — H)c colors per type. 



inconsistent with a constant slope for all spirals, given the 
uncertainties, a similar dip is also obvious in the slopes 
of the {B — V)t- and {U — -B)T-logio(^i25) relations, sug- 
gesting that this is a real phenomenon. One should first 
remember that the slope is not a measure of the total ex- 
tinction but of the reddening relative to face- on. With a 
simple radiative transfer model where the stars and the 
dust are distributed homogeneously in an infinite slab, we 
find that when the face-on optical depth ry in the V- 
band increases, the slope begins first to steepen (Fig. |^) 
but when ry becomes larger than 1-2, the slope fiattens 
although the colors go on reddening (Fig. ||a). Although 
this modeling is simplistic and we can not infer from it 
the value of the optical depth, this suggests that the ex- 
tinction may indeed be higher for Sbc and that the slope 



has begun to decrease. This conclusion is reinforced by 
the fact that the slope decreases when we consider effec- 
tive B-H colors (/3i[logio(-R25)] = 0.60/0.55), which are 
probably more extinguished than total colors. 

We observe a positive slope for irregular galaxies in 
B — H, contrary to B — V and U — B where none is de- 
tected. However, the uncertainties are large and the null 
slope is within 2a. Moreover, our "Im" type regroups in 
fact not only Im galaxies (T — 10) but also Sm (T = 9) 
which are more elongated and redder. The small size of 
the "Im" sample makes it impossible to decide whether 
this reddening is due to the dust or to intrinsically redder 
populations in Sm galaxies. 

No slope is observed for lenticular galaxies, which may 
indicate that no significant amount of dust is present in SO 
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type 


N 


M{B - 


K)t] 




f3o[{B-KU 


a 


E 


142 


3.98 ± 


0.02 





13 ± 


01 


4.08 ±0.02 





15 


± 





01 






3.95 ± 







09 ±0 


04 


4. 06 ± 0.02 







± 





02 


SO 


159 


3.97 ± 


0.02 





18 ± 


02 


4.07 ±0.02 





19 


± 





01 






3.95 ± 


0.02 





15 ±0 


03 


4.06±0.02 





19 


± 





oa 


Sa 


99 


3.79 ± 


0.04 





37 ±0 


03 


3.99 ±0.04 





36 


± 





03 






3.79 ± 







33 ±0 


04 


4.00 ± 0.04 





55 


± 





03 


Sb 


100 


3.64 ± 


0.05 





41 ± 


04 


3.84 ±0.05 





42 


± 





04 






3.62 ± 


0.05 





42 ±0 


05 


3.82 ± 0.05 





^5 


± 





0^ 


Sbc 


48 


3.47 ± 


0.08 





46 ± 


06 


3.64 ± 0.07 





40 


± 





05 






3.45 ± 


0.07 





32 ±0 


10 


3.64 ± 0.06 





32 


± 





06 


Sc 


46 


3.26 ± 


0.09 





49 ± 


07 


3.46 ±0.09 





50 


± 





06 






3.22 ± 


0.0,9 





39 ±0 


iO 


5.^5 ± 0.08 





^5 


± 





07 


Sd 


24 


2.89 ± 


0.10 





39 ± 


08 


3.09 ±0.10 





43 


± 





08 






2.80 ± 


0.08 





21 ±0 


i3 


3.01 ± 0.09 





31 


± 





JO 


Im 


25 


2.23 ± 


0.11 





49 ± 


09 


2.41 ±0.11 





50 


± 





08 






g.^i ± 


O.ii 





45 ±0 


10 


2.40 ± O.ii 





50 


± 





08 



Table 6. {B — K)t and {B — K)e colors per type. 



type 


TV 




H)t] 


/3l[logio(i?25)] 


Ml 


a 


E 


138 


3.77 ± 


0.02 


-0.77 ±0.19 


0.10 





12 


± 





01 






5.75 ± 


0.02 


-1.00 ± 0.27 




(0 


05 


± 





05) 


SO 


153 


3.72 ± 


0.02 


0.04 ±0 


11 


0.20 





17 


± 





02 






3.69 ± 


0.02 


0.04 ± 


12 







i5 


± 





05 


Sa 


110 


3.53 ± 


0.04 


1.00 ± 


22 


0.19 





31 


± 





03 






5.55 ± 


0.04 


0.87 ± 


J9 







24 


± 





0^ 


Sb 


171 


3.48 ± 


0.03 


1.03 ± 


14 


0.27 





34 


± 





02 






3.45 ± 


0.03 


1.12 ±0 


14 







54 


± 





05 


Sbc 


110 


3.29 ± 


0.04 


0.73 ± 


19 


0.25 





34 


± 





03 






3.27 ± 


0.03 


0.7a ± 


2.8 







2,9 


± 





0^ 


Sc 


125 


3.14 ± 


0.04 


0.89 ± 


15 


0.29 





37 


± 





03 






3.12 ± 


0.04 


0.95 ± 


J7 







32 


± 





05 


Sd 


135 


2.89 ± 


0.04 


0.76 ± 


15 


0.34 





40 


± 





03 






2.88 ± 


0.04 


0.82 ± 


16 







5,9 


± 





05 


Im 


41 


2.41 ± 


0.08 


0.72 ± 


36 


0.29 





40 


± 





06 






2.39 ± 


0.07 


0.70 ± 


5^ 







55 


± 





0« 



Table 7. (i? — iJ)T color per type as a function of Xn — Iog2o(-R25)- 



(see also Sandage & Visvanatan 197^ ) . Another possibility 
is that the disk is overwhelmed by the bulge. No inclina- 
tion dependence of the colors would then be expected if 
the bulge is spherically symmetric or if it is devoid of dust. 

The most striking result is the negative slope observed 
for ellipticals. This certainly does not come from the dust 
because the slope would be positive. For ellipticals, R25 is 
not directly related to the inclination because their true 
shape is unknown. In the mean, rounder ellipticals seem 
redder and are probably more metal-rich (Terlevich et al. 



198C ). They also tend to be brighter (Tremblay & Merritt 
1996) which may link the color-log]^Q(i?25) relation of ellip- 



ticals to their color-magnitude relation. No such relation 
is observed in B — V , which is anyway a poor indicator 
of the metallicity. The [U — -B)T-logio(i?25) is not very 
conclusive: whereas the ML estimator provides a negative 



slope, as in i? — , the MCES estimator does not detect 
any. Since no obvious bias between the two estimators is 
observed for the other types and the intrinsic scatter is 
very small compared to the uncertainties on the colors for 
ellipticals - which is the worst case for the MCES esti- 
mators -, we tend to trust more the result given by the 
ML. 

5.5. The color-magnitude relation 
5.5.1. Global relation 

The color-magnitude (CM) relation is an important con- 
straint on the dynamical models of galaxy formation for 
the absolute magnitude may be related to the mass of 
the galaxy via a mass-to-luminosity ratio. An important 
question is the choice of the reference band, optical or 
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Fig. 7. Global {B — NIR)T-absolute B magnitude relations. 
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Fig. 8. Global {B — NIR)T-absolute NIR magnitude relations. 
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- NIR)t 
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2.68 ±0.02 


0.03 ±0.02 
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0.47 ±0.02 


2.68 ±0.02 


-0.20 ±0.02 


-23.34 


0.41 ± 0.01 
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-0.05 ± 0.03 




0.46 ± 0.02 


2.67 ±0.02 


-0.19 ± 0.02 




0.40 ± 0.02 


(B 


-H)t 
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3.41 ±0.02 


-0.05 ±0.01 


-20.46 


0.45 ±0.01 


3.42 ±0.02 


-0.19 ±0.01 


-23.90 


0.40 ±0.01 








3.39 ± 0.02 


-0.06 ±0.02 




0.45 ± 0.01 


3.40 ± 0.01 


-0.18 ± 0.01 




0.40 ± 0.01 


{B 


-K)t 


546 


3.68 ±0.02 


-0.08 ±0.02 


-20.67 


0.49 ±0.02 


3.69 ±0.02 


0.24 ±0.01 


-24.38 


0.41 ±0.01 








3.66 ± 0.02 


-0.10 ± 0.03 




0.49 ± 0.02 


3.68 ± 0.02 


-0.23 ± 0.02 




0.40 ± 0.01 



Table 8. {B — NIR)t vs. Babs or NIRabs color-magnitude relations, where NIR — J, H or K. 



^ % ^ ^ I 



i 



Sb Sbc 
type 



-m 

E SO So Sb Sbc Sc Sd Im 

type 



* 



Sb Sbc 
type 



Fig. 5. Slopes of the {B - H)t-, [B - V)^-, and {U - 
B)^-\ogif){R25) relations as a function of type. The black 
and white squares correspond respectively to the ML and 
MCES estimators. 



NIR, used for the absolute magnitude. Usually, the NIR 
has been adopted because it suffers little extinction. An- 
other reason is that it is produced mainly by old giants 
and is therefore expected to be a better indicator of the 
mass than the optical, which is more sensitive to the re- 
cent star formation. We have plotted on Figs. and || the 
color-absolute magnitude relation with either the optical 



0.2 0,4 0,6 0.8 



T =0.3 

V 

V ^ 

T^ = 2 

T =3 

V 

T =5 



Fig. 6. {B — iJ)-log]^Q(i?25) relations computed with the 
radiative transfer code (see text) for different face-on V- 
band optical depths. We assume go = 0.18, a typical value 
for an Sb. Top: reddening relative to the dust-free case. 
Bottom: reddening relative to face-on. 



or the NIR as reference band. The corresponding estima- 
tors are given in Table |. The (B-NIR)T-NIRabs relation 
is much steeper than the {B — NIR)T-Sabs one. If there 
was no intrinsic scatter, the slopes should be nearly the 
same. For example, if (B — NIR)t — a NIRabs ±k, then we 
should have {B - NIR)t = [a/(l ± a)](Sabs + «), which 
has almost the same slope since |q;| << 1. The scatter is 
moreover smaller when the NIR is the reference band. To 
detect a color-magnitude relation, the choice of the NIR as 
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reference band is thus clearly favored. Though the overall 
agreement is very good, we note that the ML estimator of 
the slope is slightly biased when compared to the a priori 
unbiased MCES estimator. Its slope is "more positive" in 
the {B — NIR)T-i?abs relation because of the positive cor- 
relation of the errors in {B ~ NIR)t and Sabs, and "more 
negative" in the {B — NIR)T-NIRabs relation because of 
the negative correlation of the errors in {B — NIR)t and 

NIRabs. 

The slope of the {B — NIR)T-absolute magnitude re- 
lation increases with the wavelength of the NIR band for 
there is also a CM relation in the NIR colors. This effect 
is less apparent in [B — H)t because the B — H sam- 
ple contains many late-type galaxies which have a flatter 
slope: a linear fit for all the galaxies together is hence too 
simplistic. 

5.5.2. Relation per type 

The global color-magnitude relation discussed above is ac- 
tually dominated by star-forming galaxies (/3i(i?abs) = 
-0.21/ -0.20 against -Q.Q7/-0.05 forE/SO). These val- 
ues are close to the slopes of the — K)[D25]-Ki^\^s re- 
lations (-0.24 ± 0.03 for Sa-Sdm and -0.09 ± 0.04 for 



E/SO) determined by Mobasher et al. ( |1986[ ). The slope 
for spiral galaxies is much steeper than the relations ob- 
tained by Bershady (1995) for its spectral types bk, bm, 
am and fm. This is however not surprising, because using 
spectral types rather than photometric types bins galaxies 
as a function of their colors and tends to suppress any 
color-magnitude relation. 

An especially interesting question is whether the mass 
is the main driving force of galaxy evolution and star for- 
mation, or if the galaxy type has also to be considered. 
To answer this question, it is worth to look at the color- 
magnitude relation per type more thoroughly. We have 
plotted the {B — 11)^-11^x^5 relation as a function of type 
in Fig. |. We also show in Fig. ^ the [B - H)^-H^hs 
relations where {B — H)^, is the color corrected to face- 



on using the /3i[log2o(-R25)] determined in |5.4 We assume 
that the extinction is negligible in the NIR and do not cor- 
rect ifabs- No correction is also applied to elliptical colors 
because their color- i?25 relation is not due to the extinc- 
tion. The values of the estimators are given in Table ^. 
We note that the intrinsic scatter in the corrected colors 
is significantly reduced. A more picturesque comparison of 
the relations for the different types is plotted on Fig. |ri|. 
Though the agreement between the estimators is not as 
good as previously and the uncertainties are large, the 
slope obviously depends on the type. We tentatively sum- 
marize the following results: 

Star-forming galaxies: the slope is much fl atter for Sd-Im 
than for Sa-Sc galaxies. Mobasher et al. ( 1986 ) probably 
did not detect this because the few Sd/Im galaxies in their 



late-spirals sample were lost among Sbc and Sc galaxies. 
At a given NIR absolute magnitude, the mean colors of the 
types are different, indicating that the CM relation is not 
simply a type-mass relation: the NIR intrinsic luminosity 
must be completed by the type to characterize the colors 
(and the star formation history) of star- forming galaxies. 
A similar conclusion was drawn by Gavazz i ( |1993D but was 
later challenged by Gavazzi et al. ( |1996b[ ). 

E/SO: the color-magnitude relation of early-type galax- 
ies is usually explained by the galactic winds produced 
by super- novae in a starbursting environment, which ex- 
pel th e gas and quench the star formation (Matthews & 
Baker 1971 ). Massive (and NIR bright) galaxies having a 
deeper gravitational well, they are able to retain the gas 
longer and to prolong the star-forming phase. This results 
in a higher mean stellar metallicity and thus in redder 
colors (Faber 1973 ). Besides its interest to constrain the 
models of galaxy formation, the CM relation has also been 
proposed as a distance indicator. 

We find a different slope for ellipticals and lenticulars. 
The slope of SO seems closer to that of early and inter- 
mediate spirals than to that of ellipticals, which may give 
some clues on their evolution. 

The MCES estimation does not detect any significant 
slope for ellipticals. Because of the small intrinsic scatter, 
it is maybe not the best estimator^, but even the slope 
determined by ML is smaller than that of the compara- 
ble {V — K) [5/i~^ kpcJ-Vr relation established by Bower et 
al. ( |1992| ) in the Virgo and Coma clusters. Our intrinsic 
scatter is also larger that theirs. A possible explanation 
of this discrepancy is that our sample contains not only 
cluster galaxies but also field ellipticals, which may have 
a different star form ation history and increa se the scatter 
(Lars on et al. 198C; Kauffmann & Chari ot 1998 ; Baugh 
et al. 1996| ; see yet Bernardi et al. 199S and Schade et 
al. 1996). However, we also note that the relation estab- 



lished by Bower et al. is based on aperture colors. Their 
aperture is comparable to the mean effective aperture of 
ellipticals in Virgo and samples only about half their lu- 
minosity (cf. also Kodama et al. 1997). Since the effective 
aperture of bright ellipticals is larger, only their inner re- 
gions are observed in the Bower et al.'s aperture. Because 
of the blue-outwards color gradient, their aperture color is 
redder and the slope is higher than the one derived from 
total colors. We therefore encourage the modellists of el- 
liptical galaxies, either to compare their relations to total 
observed colors, or to predict aperture colors. 

The color-magnitude relations as a function of the B 
absolute magnitude are finally plotted for each type on 



^ Note that the slopes obtained for the effective color 
vs. absolute magnitude relation with the two estimators 
(—0.03/— 0.0^) are in good agreement, both between them 
and with the MCES slope of the total color-magnitude rela- 
tion. 
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Fig. 9. {B — iJ)T-absolute H magnitude relations per type. 
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Fig. 12. {B — i/)T-absolute B magnitude relations per type. 
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type 


N 


f3ol{B-H)T] 


/3l(-ffabs) 


Ml 










E 


119 


3 


79 ± 0.02 


-0.05 ± 0.01 


-24.95 


0.12 ± 0.02 




3.79 ±0.02 


0.05 ± 0.01 


0.12 ± 0.02 






3 


76 ± 0.02 


-0.02 ± 0.02 




{0.06 ± 0.0^) 




3.76 ± 0.02 


-0.02 ± 0.02 


{0.06 ± 0.0^) 


SO 


110 


3 


70 ± 0.02 


0.13 ± 0.02 


-24.00 


0.14 ±0.02 




3.67 ±0.02 


0.12 ± 0.02 


0.18 ± 0.02 






3 


67 ± 0.02 


-0.10 ± 0.05 




0.14 ± 0.05 




3.64 ± 0.02 


-0.11 ± 0.03 


0.16 ± 0.03 


Sa 


94 


3 


52 ± 0.04 


0.22 ± 0.04 


-23.85 


0.37 ±0.03 




3.35 ±0.04 


0.20 ±0.03 


0.29 ± 0.03 






3 


51 ± 0.0^ 


-0.22 ± 0.05 




0.36 ± 0.0^ 




3.35 ± 0.0^ 


-0.20 ± 0.05 


0.28 ± 0.0^ 


Sb 


140 


3 


53 ± 0.04 


0.14 ±0.04 


-24.31 


0.37 ±0.03 




3.18 ±0.03 


-0.17 ± 0.03 


0.30 ± 0.02 






3 


51 ± 0.0^ 


-0.12 ± 0.05 




0.36 ± 0.03 




3.17 ±0.03 


-0.14 ± 0.0^ 


0.29 ± 0.03 


Sbc 


91 


3 


36 ± 0.04 


0.19 ± 0.04 


-24.12 


0.35 ± 0.03 




3.14 ± 0.03 


—0.24 ± 0.03 


0.26 ± 0.03 






5 


33 ± 0.0^ 


-0.12 ± 0.07 




0.32 ± 0.0^ 




3.i2 ± 0.03 


-0.18 ± 0.06 


0.25 ± 0.03 


Sc 


102 


3 


22 ± 0.04 


0.19 ± 0.04 


-23.44 


0.39 ± 0.03 




2.87 ±0.03 


0.23 ±0.03 


0.27 ±0.03 






3 


19 ± 0.04 


-0.17 ± 0.05 




0.38 ± 0.03 




2.85 ± 0.03 


-0.21 ± 0.0^ 


0.25 ± 0.03 


Sd 


95 


2 


98 ± 0.05 


-0.08 ± 0.04 


-22.11 


0.46 ± 0.04 




2.65 ±0.04 


-0.13 ±0.03 


0.37 ±0.03 






2 


95 ± 0.05 


-0.07 ± 0.05 




0.46 ± 0.0^ 




2.64 ± 0.0^ 


-0.J2 ± 0.04 


0.3^ ± 0.03 


Im 


20 


2 


55 ±0.09 


-0.07 ±0.05 


-21.60 


0.34 ±0.07 




2.32 ±0.08 


0.09 ±0.04 


0.27 ±0.06 






2 


54 ± 0.08 


-0.06 ± 0.07 




0.31 ± 0.09 




2.32 ± 0.08 


-0.08 ± 0.06 


0.24 ± O.ii 



Table 9. {B — H)t and {B — H)j, per type as a function of Xn = i/abs- 



type 


N 


/3o[{B-H)t] 


/3l(J/abs) 


Ml 


/32[logio(i?25)] 




a 


E 


115 


3.80 ±0.02 


0.05 ±0.01 


-24.97 


-0.67 ±0.19 


0.10 


0.10 ± 0.02 






3.76 ± 0.02 


-0.02 ± 0.02 




-0.S3 ± 0.35 




{0.04 ± 0.03) 


SO 


106 


3.66 ±0.02 


-0.16 ±0.02 


-24.05 


0.37 ±0.16 


0.19 


0.18 ± 0.02 






3.66 ± 0.02 


-0.i2 ± 0.04 




0.27± 0.^7 




0.14 ± 0.03 


Sa 


94 


3.47 ±0.04 


0.23 ±0.03 


-23.85 


1.18 ± 0.23 


0.17 


0.30 ± 0.03 






3.49 ± 0.04 


-0.21 ± 0.04 




1.21 ± 0.20 




0.27 ± 0.04 


Sb 


140 


3.47 ±0.03 


0.17 ±0.03 


-24.31 


0.97 ±0.13 


0.25 


0.30 ± 0.02 






3.45 ± 0.03 


-0.14 ± 0.0^ 




0.99 ± 0.i3 




0.29 ± 0.02 


Sbc 


91 


3.30 ±0.03 


-0.28 ±0.03 


-24.12 


1.22 ± 0.16 


0.25 


0.23 ± 0.03 






3.28 ± 0.03 


-0.22 ± 0.06 




1.19 ± 0.2i 




0.24 ± 0.03 


Sc 


102 


3.13 ±0.04 


-0.24 ±0.03 


-23.44 


1.10 ± 0.13 


0.29 


0.27 ±0.03 






3.Ji ± 0.03 


-0.22 ± 0.0^ 




1.17 ± 0.14 




0.24 ± 0.03 


Sd 


95 


2.91 ±0.04 


0.14 ±0.03 


-22.11 


1.15 ± 0.16 


0.32 


0.33 ±0.03 






2.88 ± 0.0^ 


-0._Z5 ± 0.04 




1.23 ± 0.i6 




0.33 ± 0.03 


Im 


20 


2.52 ±0.07 


0.07 ±0.04 


-21.60 


1.11 ± 0.29 


0.25 


0.18 ± 0.05 






2.53 ± 0.06 


-0.05 ± 0.05 




1.09 ± 0.29 




{0.04 ± 0.06) 



Table 10. {B — H)t color per type as a function of Xn = i?abs and Xi2 = logiQ(i?25)- 



figure The slopes are very different of those obtained 
when the NIR is the reference band. For spiral galaxies 
especially, they tend to have a positive sign, with brighter 
galaxies (in B) being bluer! This behavior is typically ex- 
pected if the relation is dominated by the extinction. If we 
assume that the extinction is negligible in the NIR, the 
same correction applied to colors may be used to com- 
pute face-on absolute magnitudes i?abs- ^^'^ {B — H)^,- 
B^^^ relations per type show then more reasonable slopes 
(Fig-H). 



5.6. The color vs. R25 & absolute magnitude relation 

We finally give in Table |l^ the estimators for a simul- 
taneous fitting of the {B — H)t colors as a function of 
log]^Q(i?25) and -ffabs in each type. The overall agreement 



with the slopes determined previously from the separate 
fitting of the colors with respect to either R25 or iJabs 
is correct. The requirement that ffabs be known tends to 
select the brightest galaxies and may explain the small 
discrepancies in the slope clS 0(S Qi function of logjQ(i?25), 
as this one depends probably also on the mass {^ ffabs) of 
the galaxy. The very small scatter obtained for Im galaxies 
indicates that there are too few of them for the simulta- 
neous analysis. 

6. Discussion 

NIR growth curves of the magnitude as a function of the 
aperture have been built from the Catalog of Infrared Ob- 



servations (Gezari et al. 1997) and optical catalogs. Using 



these, we have been able to compute total NIR apparent 
and absolute magnitudes, NIR and optical-to-NIR colors. 
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Fig. 11. Analytical {B ~ i?)!p-absolute iJ magnitude re- 
lations. The segments correspond to the mean of the 
ML and MCES fits and extend from /io(^?abs) + 1-5 to 
A*o(i?abs) - 1-5 for each type. 

and to estimate their uncertainties for a large sample. A 
statistical analysis of the colors as a function of type, in- 
clination and luminosity, using estimators taking into ac- 
count the uncertainties in the variables and the intrinsic 
scatter in the colors, highlights the interest of the NIR and 
notably of the comparison of the optical to the NIR. 

Optical-to-NIR colors show a well defined sequence 
with type: the mean irregular is 1.3 magnitude bluer in 
B — H than the mean elliptical. The intrinsic scatter is 
higher for star-forming galaxies {a ^ 0.4) than for ellipti- 
cals or lenticulars {a ~ 0.1-0.2). 

Because of the small extinction in the NIR, the optical- 
to-NIR colors of spiral galaxies redden considerably with 
increasing inclination, putting therefore constraints on the 
amount of the dust and the respective distributions of 
dust and stars. SO colors do not depend on the inclina- 
tion. Rounder elliptical are redder than more elongated 
ones. 

The color-absolute magnitude relation is much steeper 
and tighter when the NIR is used as the reference band 
rather than the optical. Examination of the color vs. B 
absolute magnitude relation for each type suggests that 
this is due to the extinction which dims and reddens 
the galaxies and counteracts the color-NIR relation. A 
color-magnitude relation exists in each type, with brighter 
galaxies in the NIR being redder. The slope is steeper for 
Sa-Sc than for early-type galaxies or Sd-Im. The relation 
we obtain for ellipticals is shallower than in other stud- 
ies. A likely explanation is that this is due to the effect 
of the color gradient in small apertures, a lesser, i.e. red- 
der, fraction of the galaxy being observed in large bright 
ellipticals than in small faint ones - a problem we do not 
have using our total (or we expect so) colors. The slope 
for SO galaxies is closer to that of spirals than to that 
of ellipticals. The color at a given absolute magnitude be- 



comes bluer with increasing type, indicating that both the 
mass and the type must be used to characterize the col- 
ors and to describe the star formation history of galaxies. 
Once corrected for the inclination and the magnitude ef- 
fect, the intrinsic scatter in the colors of spiral types drops 
significantly, although a few outliers with very blue colors 
typical of starbursting galaxies remain. 

However satisiying these results are, some problems 
still exist. NIR observations are obtained in small aper- 
tures and the extrapolation to total magnitudes and col- 
ors suffers from significant uncertainties. Nevertheless, we 
do not believe that our results strongly depend on this ex- 
trapolation since similar slopes are obtained when we use 
effective colors - which are interpolated or only slightly 
extrapolated - instead of total colors. 

The selection of galaxies is a more crucial problem. 
The only criterion used here has been to take whatever 
was available in the NIR. For this reason, the sample, 
especially in J and contains few late-type galaxies 
but many cluster galaxies, Seyferts and other active nu- 
clei. We also lack of intrinsically faint galaxies. The color- 
magnitude relations we have determined should therefore 
not be extrapolated at i?abs fainter than ~ —20. Large 
samples as the Sloan Digital Sky Survey in the optical or 
DENIS and 2MASS in the NIR should solve this problem. 

The last important problem is the bias between the 
two estimators we use, especially in the slope of the color- 
magnitude relation. The discrepancy tends to disappear 
however when we consider effective colors which have 
smaller uncertainties. The MCES estimators seem to be 
closer to the true slopes than the ML, but the latter pro- 
vides a better estimate of the intrinsic scatter and has a 
lower variance. 

Combining the NIR and the optical enlarges consid- 
erably our vision of the stellar populations and the dust 
content of galaxies. In our next paper, we will extend our 
analysis to the near ultraviolet-optical wavelength range. 
We will publish "color" energy distributions of galaxies 
from the near-UV to the NIR that will be used as tem- 
plates to constrain the star formation history of galaxies 
of different types and masses. A further step will be the 
extension of these techniques to the far-UV and the far-IR. 

Acknowledgements. This research has made use of: 

— the Catalog of Infrared Observations operated at th e 
NASA/Goddard Space Flight Center (Gezari et al. |l997| ), 

— the Hypercat catalog operated a t the CRAL-Observatoire 
de Lyon (Prugniel & Heraudeau 1998), 

— the Third Refere nce C atalog of Bright Galaxies (de Vau- 



couleurs G. et al. 1991) 



the NASA/IPAC Extragalactic Database (NED) which is 
operated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. 

We also acknowledge: 



M. Fioc & B. Rocca-Volmerange: Optical-to-NIR colors of nearby galaxies 



23 



the Statistical Consulting Center for Astronomy operated 
at the Department of Statistics, Penn State University, 
M.G. Akritas (Director), 



Numerical Recipes by Press et al. (1992) 



M. F. acknowledges support from the National Research 
Council through the Resident Research Associateship Pro- 
gram. 



Appendix A: The optical growth curve 

Let Sn{X) be defined by 



Sn{X) = -2.5 log 



10 



f 

Jo 



exp(— a;) dx 



Jq°° x^" 1 exp(— x) dx 



where Xq is such as iS„(0) = 2.51ogio(2) = 0.7 5. 

According to Prugniel & Heraudeau ( |1998| ), B{X,Tp 
is defined by: 

BiX,Tp) - Sn{X), n = (15 - Tp)/5 



if Tp < —5 or Tp > 10 - corresponding to a Sersic ( |1968D 
luminosity profile i{A) oc exp [— (yl/Ao)^/"] and by 

B{X,Tp) = (2/3-rp/15)54(X) + (l/3 + rp/15)5i(X) 

if — 5 < Tp < 10, i.e., interpolated between the de Vau- 
couleurs profile and the exponential profile valid respec- 
tively for giant ellipticals and pure-disk galaxies. 



Appendix B: Computation of NIR magnitudes 
and uncertainties 

Let us first define Bi = B{Xi, Tp), 

n 

S ^^W^, 
i=l 

where Wi is the weight attributed to the i**^ point, 

n 

Sb = y^^WjBi, 

i=l 
n 

Sb^ = ^w,Bf, 

4=1 

n 

S,n = y^^WjUli, 
i=l 

n 

SBm = WiBiUli, 

1=1 

and 

A = SSb^ — 5*0. 

Two methods have been used, depending essentially on 
the number of observations and their accuracy. In the first 



one, the total NIR magnitude mx is computed assuming 
s = So by minimizing 



Wi (rrii - mx - soBi)"^ 

i=l 

with respect to mx- We obtain 

mJp'' {Srn - SoSb)/S, 



,(1) 



i=l 



and the corresponding effective magnitude is 
=mx + 2.51ogio(2)so, 



< [2.51ogio(2)5 - Sb? +Y.w^a^ /S\ 



We assume here that, for want of constraint on s, the 
uncertainty on s (o'(s)) is equal to the intrinsic scatter 
a so for this type. 

In the second method, s is considered as a free param- 
eter and mx is computed by minimizing 

n 

Wi {rrii - mx - sBi)^ 

i=l 

with respect to mx and s. We then get 
m^'' = (SB^Sm - S'eS'B„i)/A, 



,(2) 



.i=l 



/A^ 



S^'^ = (SSBrn - 5b5„,)/A, 



.(2) 



,(2) - 



mx + 2.51ogio(2)s, 



/A^ 



and 



,(2) 



. i=l 



^wfcrf [Sb2 - B^SB 

+ 2.51ogio(2) iSB,-SB)f\ /A^ 



Note that we have made the questionable assumption that 
the errors in the aperture magnitudes of a given galaxy are 
independent. 

We are however not primarily interested in the best 
fitting growth curve, which would be determined assum- 
ing Wi = l/cf, but rather in the best estimate of the 
asymptotic magnitude. When there is only one point, the 
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uncertainty on mx is a{m,T) = {(jf + Sicr(s)^)^/^, i.e., the 
uncertainty on tot decreases with the size of the aperture. 
It is therefore reasonable to adopt Wi — l/{crf + Bfa{s)'^) 
for all the points to give more weight to the large aper- 
tures. 

To deal with outliers, we apply an iterative proce- 
dure. We initially assume ct(s) = asg for the slope and 
Ci = (Jmo for all the points. At each step, we fit the 
growth curve to the data, compute a{s) from above equa- 
tions, and estimate new values of the uncertainties by 
(Tj = max((Tmo: I'^i —mT — sBil). This procedure is a com- 
promise between our a priori uncertainty ^'Hd the a 
posteriori estimate Er=i("^« ~ sBiY /{n — fc)]^/^, where 
A: = 1 or 2 is the number of fitted parameters (and also 
the index of the method!). It has the advantage of reduc- 
ing the weight of outliers and takes the scatter around 
the growth curve automatically into account. Practically, 
it improves the fit and gives a reasonable estimate of the 
uncertainties. The convergence is usually achieved in a few 
iterations. 

The second method provides a better fit than the first 
one but is less secure because the slope of the growth 
curve is free. It is used only when following conditions are 
fulfilled simultaneously: 

1. n > 3, which allows to check the validity of the fit. 



2. a 



.(2) 



3. and 



(1) 



" r 

< n^-f 



i.e., the fit provided by (2) is better enough than (1) 
to justify the deviation of s relatively to Sq. 
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